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ABSTRACT 


We present detailed calculations of accreting magnetized neutron star 
atmospheres heated by the gradual deceleration of protons via Coulomb 
collisions. Self-consistent determinations of the temperature and density 
structure for different accretion rates are made by assuming hydrostatic 
equilibrium and energy balance, coupled with radiative transfer. The full 
radiative transfer In two polarizations, using magnetic cross sections but 
with cyclotron resonance effects treated approximately, is carried out in the 
Inhomogeneous atmospheres. For ft < 10 17 g s” 1 , the equilibrium atmospheres 
have temperatures and optical depths which are very sensitive to the strength 
of the surface magnetic field. Because of a decreased efficiency of cyclotron 
line cooling, atmospheres with higher magnetic fields are hotter, more 
optically thin, and radiate harder spectra. The computed pencil beam pulse 
shapes show frequency dependent structure: for a wide range of aspect angles, 

the pulses switch from single to multiple below - w H /4, and there is a general 
hardening of the spectral index toward midpulse. 
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I. INTRODUCTION 

An outstanding problem in the theory of accreting magnetized neutron 
stars is the lack of a unique, self-consistent model of the accretion column, 
resulting In major uncertainties about how the spectrum and the pulses are 
produced. Given that many of our determinations of neutron star masses, radii 
and magnetic field strengths have come from such accreting systems, it is 
vital to understand in more detail the origin of the radiation from which all 
this information is gleaned. It is known that for high luminosities, 

L > 10 37 erg s— 1 , radiation pressure plays a major role in decelerating the 

A 

accreting matter (Basko and Sunyaev 1976, Wang and Frank 1981). However for 
luminosities L < 10 37 erg s — x , to which we shall confine ourselves, it is not 

A 

clear whether the deceleration occurs gradually via Coulomb encounters, 
through inelastic nuclear collisions, or by the formation of a collisionless 
shock (Zeldovich and Shakura 1969, Lamb, Pethick and Pines 1973, Langer and 
Rappaport 1982). This question is still unsettled, and because of this basic 
uncertainty the best approach at present is to investigate the observational 
consequences predicted by different models, in the hope of thus being able to 
distinguish among them. 

In the present paper we study the gradual deceleration model, and argue 
that in most cases Coulomb collisions will be more effective than nuclear 
collisions. We calculate the temperature and density gradients in hydrostatic 
and energy equilibrium, coupled with the fully macnetized radiative transfer 
including vacuum polarization effects,- Preliminary results of this 
calculation have been presented in Meszaros et al. (1983). We include the 
effects of cyclotron line cooling in an approximate way, but do not attempt to 
discuss the details of the line spectrum, since a full theory for the transfer 
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near the resonance Is still missing (Wasserman and Salpeter 1980, Bussard 

' t 

1980, Kirk and Meszaros 1980, Nagel 1980). We find that such Coulomb heated I 

i 

atmospheres are rather shallow, of height z < 10 3 cm which is much less than 

the radius of the emitting polar cap. Previous calculations of the spectrum 

for such plane parallel atmospheres with $ parallel to z have previously been 

made assuming a homogeneous slab or semi -infinite medium, with an estimated j 

density and temperature (Bonazzola et al 1979, Ventura et al 1979, Meszaros et 

.1 

Jj 

al 1980, Nagel 1980). Here we calculate the continuum spectrum from the self- 
consistent density and temperature stratification of the medium. We have not f 

i* 

Included the possible effect of absorption by matter farther away than the j 

immediate vicinity of the polar cap regions (McCray and Lamb 1976, Basko and 

| 

Sunyaev 1976), nor the possible effects of interstellar absorption, which may I 

ij 

cause a low energy turnover. Incoherent scattering is implicit in our 

treatment of the line transfer, and Compton cooling is included in the energy | 

> 

balance and in the atmospheric structure calculation. However, the continuum 

spectrum is calculated assuming coherent scattering. Introducing incoherent \ 

tr 

i 

scattering should not c^nge the continuum spectrum very much for our range | 

f 

of ft, but an accurate prediction of observable cyclotron line profiles would j 

require a more elaborate treatment. 

I 

We have also calculated the intrinsic beaming pattern of the radiation 4 • 

escaping from the stratified atmosphere. Previous calculations of beam f j 

patterns for plane parallel atmospheres with § perpendicular to the surface 
assumed an estimated homogeneous density and temperature (Tsuruta 1974, Basko 
and Sunyaev 1975, Nagel 1981a, Meszaros and Bonazzola 1981). Calculations of 
beam patterns from atmospheres with cylindrical geometry (Yahel 1980, Nagel 
1981a, Pravdo and Bussard 1981) have also assumed a homogeneous density and 
temperature, or taken an estimated distribution. The beam patterns found in 
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the present work, arising from the self-consistent inhomogeneous (plane 
parallel) atmospheres heated by Coulomb deceleration are of the pencil beam 
type. They are characterized at many aspect angles by single pulses at high 
frequencies which break up into multiple pulses at lower frequencies. The 
breakup frequency, within some limitations, corresponds to a submultiple of 
the cyclotron frequency, and this can give an Indirect estimate of the field 
strength. We also discuss the phase dependence of the spectrum, in terms of 
spectral index vs. phase plots. The present models show a general, but not 
universal tendency for the spectrum to change most rapidly toward midpulse. 

The details of both the pulse shape and the spectral index vs. phase behavior 
are dependent on the magnetic field strength and the viewing angle, as well as 
the specific temperature and density gradients in the atmosphere. 

In the final section of the paper we discuss some of the implications of 
these calculations, and compare them with other models and with available 
observations. 

II. PHYSICAL PROCESSES IN STRONG MAGNETIC FIELDS 

There are several major complications which arise in the presence of a 
strong magnetic field due to the quantization of electron energy perpendicular 
to the field. First, the radiation propagates in tv/o normal polarization 
modes, each with rather different opacities, which are coupled by polarization 
exchange scattering. Second, the proton stopping is somewhat more 
complicated, as the strong field will partially inhibit the exchange of 
momentum with the atmospheric electrons. The equations are of a similar 
character to the nonmagnetic case as treated by Zeldovich and Shakura (1969) 
and Alme and Wilson (1973), but a detailed treatment of the magnetic 
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ORIGINAL PAtte 
OF POOR QUALITY 

peculiarities makes it necessary to resort to a numerical calculation. 

We assume a plane parallel atmosphere, varying only along the z 
coordinate perpendicular to the surface (i.e. parallel to the field $). The 
grammage traversed by a proton falling into the atmosphere is defined by 

y = p (z) dz (1) 

in units of g cm"^. The plane parallel atmosphere is assumed to extend 
sideways over an area A - tt R 2 , where R is the magnetic polar cap radius, 

r r 

«• 

there being one such area at each pole. For a cl. lain magnetic field strength 
and accretion rate, we determine the polar cap radius R p as a function of the 
Alfven radius. For an Alfven surface interacting with an accretion disk, we 
use a formula originally due to Lamb et al (1973) for the Alfven radius, which 
we rewrite as 

18/69 -16/69 -13/69 120/69 40/69 

R A (di sk ) = 3 x 10 7 a rti m Rg B.j 2 (cm) (2) 

We have used the disk parameters of Shakura and Sunyaev (1973), case (b), 
where a T > p gas > p rad , and a <_1 is the viscosity parameter, m = M/M 0 , 
m = fa/1.3 x 10 18 m, R g - R N /10 6 cm, B 12 = B/10 12 G. In. terms of R A and R^, 
the polar cap radius is (Rp/R^) 2 - R^/Ra* so that 

-9/69 8/69 13/138 -51/138 -20/69 

(R p /R N ) - °-l 6 a itim Rg B 12 (3) 

Thus, for B = 4.4 x 10^ G, and M = 10 17 , 10 16 , 10 15 and 10 1 '♦ g r 1 , m = a 
=1, R 6 = 1.2, we get respectively R p = 8.83, 6.66, 5.0 and 4.6 x 10 4 cm. 
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a. Deceleration of the accretion flow 


The Interaction between the accreting material and the atmosphere was 
treated using the techniques described by Kirk and Galloway (1981, 1982), In 
which the effect of small -angle Coulomb scattering is considered. Early work 
in this field (Basko and Sunyaev 1975, Pavlov and Yakovlev 1976) had Indicated 
that such small -angle scattering processes were Inhibited by the strong 
magnetic field in a neutron star's atmosphere. In this case nuclear 
collisions between accreting protons and protons of the atmosphere would 
become Important, and limit the penetration depth to about 50 g cnf 2 . 

However, these calculations assumed that the accreting protons move on a 
rectilinear path through a cold electron gas. They are no longer valid when 
the thermal velocity of electrons in the atmosphere ceases to be negligible 
compared to the speed of the accreting particles, and when the gyroradius of 
the proton is of the same order of magnitude as the Debye screening length. 
Under such circumstances, the effect of small-angle Coulomb collisions between 
accreting protons and electrons of the atmosphere is, in fact, enhanced, and 
the upper limit to the penetration depth provided by nuclear collisions Is not 
always relevant. * As in the zero magnetic field case (Sitenko 1967), proton- 

*As an example, a proton of initial speed c/2 is stopped in a uniform plasma 
of n e = 10 23 cm" 3 , T e - 20 KeV, B = 5 x 10 12 gauss in a length of y 0 = 20 g 
cm’ 2 . 


proton Coulomb collisions are unimportant in the deceleration problem. 

The method of Kirk and Galloway consists of calculating the dynamic 
friction and diffusion coefficients for a heavy, char " test-proton moving in 


a strongly magnetized, homogeneous electron gas of specified temperature and 
density. The main assumptions Involved In this procedure are that the 
electron gas responds linearly to the presence of the test particle and that 
Its state is well approximated by a gas In which only the lowest Landau level 
is populated (Kirk 1980). These conditions are fulfilled In the case of 
accretion onto a neutron star provided that the density of the atmosphere 
exceeds that of the accretion flow, and that the energy density In radiation 
at the cyclotron resonance frequency Is not too large. The latter requirement 
arises because colllslonal processes are ineffective at exciting electrons 
Into higher Landau levels at densities less than 10 26 cnf^. As a result, the 
population of these levels is controlled by the radiation field, and is not 
directly related to the dispersion of electron velocity along the field 
lines. Thus, even a plasma at a temperature of 10 9 K will contain electrons 
In only the lowest Landau level, provided that the number density of photons 
at the resonant frequency Is much smaller than the corresponding blackbody 
density at a temperature T = tiw^j/k. The coefficients thus calculated are 
Inserted Into the Fokker-Planck equation, which can theft be solved to give the 
evolution in time of a specified distribution of test particles. The 
interactions of test particles with each other are neglected. The initial 
distribution is chosen to describe protons moving rapidly along the magnetic 
field lines in a homogeneous plasma. The initial distribution function assumed 
for this calculation is a Maxwellian at temperature T.* 8 500 keV centered on 
Vj _=0, ''jj ~ c / 2. Averages over the distribution at 

subsequent times then give the energy and momentum which is delivered by the 
accreting protons to the plasma of the atmosphere during the various stages of 
deceleration. 

In order to adapt this method to the inhomogeneous atmosphere of a 
neutron star, a three dimensional table was established in electron density, 
electron temperature and proton kinetic energy. At each point the rate of 
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loss of momentum dri/dy and energy dE/dy of the accretion flow was evaluated 
from one of the solutions of the Fokker-Planck equation, enabling 
interpolation over the range n e = 10*® to 10 2 ® &T >3 , T e 3 5 to 90 keV and 
inf all kinetic energy E 3 120 to 5 MeV. The dependence on B is very weak 
compared to the dependence on n e and T, since B appears only in the Coulomb 
logarithm {see Langer and Rappaport, Eqn [C4]). The neglect of the magnetic 
field strength dependence of the loss rates therefore does not introduce much 
error. Figure 1 shows the proton energy loss rates calculated in this way for 
n e 3 10 23 cnr 3 with the corresponding non-magnetie loss rates (cf. Alme amd 
Wilson 1973) also shown for comparison. 2 


2 The energy and momentum loss rates used in this paper are a factor of two 
smaller than those presented by Kirk and Galloway (1981, 1982), thus 
cancelling an error of this factor which was present in their original 
calcul ations. 


In such an interpretation it is implicitly assumed that the interaction 
between accreting electrons and atmospheric electrons proceeds relatively 
quickly. Charge neutrality is preserved by the establishment of a small 
return current in the atmosphere, in a manner similar to that in the X-ray 
producing regions associated with solar flares (Emslle 1980). Although It 
cannot be ruled out that a collisionless shock front will establish Itself 
where the accreting matter meets the atmosphere (Langer and Rappaport 1982), 
this nevertheless seems unlikely, in view of the stabilizing influence of the 
magnetic field, and the inefficiency of electrostatic instabilities (McKee 
1970, Alme and Wilson 1973, Kirk and Trumper 1982). 

In addition, we also include the effects of stopping by nuclear 


collisions, for which the energy loss rate Is taken to be 


f dE l 

WN 



(Basko and Sunyaev 1975), where the total stopping length, y N * 50 g cm' 2 . 

This deceleration rate and the corresponding momentum loss rate, taking 
p = /‘ZmF, were added to the Coulomb loss rates, since nuclear scattering Is an 
Independent process, 
b. Radiative Cooling 


We Include three contributions to the cooling of the atmosphere: 
bremsstrahlung emission Ag, Compton scattering Aq* and colli slonal excitation 
of electrons to the first excited Landau state followed by the emission of 
cyclotron line photons, A|_. We thus write the total cooling rate as 

A = Ag + Aq + A^. (4) 

The inverse bremsstrahlung (absorption) and Compton (photon energy loss to 
electrons) processes result in heating of the atmosphere and are therefore 
Included in the net cooling rate. The inverse of A L , absorption of a Tine 
photon followed by collisional de-excitation, should be a negligible 
contribution. Since the radiative de-excitation rate is much larger, line 
photons are expected to escape into the adjacent continuum by non-resonant 
scatterings before they are absorbed. 

A correct treatment of the line transfer is extremely time consuming when 
computing an inhomogeneous atmosphere by iterations, as is done below, so that 
some approximation is called for* The scattering cross section is resonant, 
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which would Increase the Inverse Compton cooling rate, but on the other hand 
photons should be scattered out of the line fairly quickly, so that most of 
the scattering should occur In the continuum. The theory of resonant line 
transfer Is extremely complicated (Wasserman and Sal peter 1980, Kirk and 
Meszaros 1980, Nagel 1981b, Bussard and Lamb 1982), and not yet fully 
understood. 

The contribution of the cyclotron resonance to the cooling is therefore 
treated In an approximate way by including a term for the production r'ate of 
line photons which is proportional to the rate of colli si onal excitations. 
Using the combined rates of electron-electron and electron-ion collisions from 
Langer and Rappaport (1982), we write the cyclotron line emissivity as 


A, « 1.22 x 1(T 20 p 


B^ 2 

esc °12 


['i.50 + 3.73 


- m.c 2 V 2 

] GXp { —j^y [ { 1+ . 04531 Bjg) "1]} 


where P esc is the probability that the line photons leave the resonance before 
being absorbed. These photons are thus assumed to appear with probability 
P esc in the line wings, where they scatter as continuum photons. As argued by 
Langer and Rappaport (1982), P esc is very near 1 for the expected densities 
and temperatures, l.e. a line photon will typically suffer a non-resonant 
scattering into the wings of the resonance before it can undergo enough 
resonant scatterings to have a significant probability of suffering a true 
absorption (radiative excitation followed by colllsional deexcltation). 
Therefore, we take P esc = 1 in almost all cases, although we have considered a 
few cases where P esc * 0.1. Actually, a value of P esc < 1 can also mimic the 
effect of a departure of the electron distribution from a pure Maxwellian. 


TV 


Langer, McCray and Baan (1980) find that departures from a Maxwellian 
distribution result In lower line emlsslvlty. The above emlsslvlty Is divided 
equally between two adjacent frequency bins, + Ao^, above and below the 
resonance frequency In both polarizations. We take the width of the line to 
be defined by the frequencies above and below the resonance where the opacity 
for resonant scattering ( Eqn [17] of Langer and Rappaport, 1982) Is 
approximately equal to the Thomson opacity. The contribution to the 
emlsslvlty from cyclotron line photons In the radiative transfer calculation 
Is therefore, 

4< , = i y ( 6 ) 

J \>L Taw^ ' ' 

where 1 = 1, 2 the two polarization modes, We assume that both modes are 

resonant and tfiat half of the cyclotron line photons are produced In each 
mode. 

Expressions for the Thomson and bremsstrahlung processes In a strong 
magnetic field have been calculated by Meszaros and Ventura (1979), Ventura 
(1979) and Nagel (1980). These are needed in the first place for finding the 
temperature and density structure of the atmosphere (radiative cooling), and 
in the second place for finding the spectrum. As In the non-magnetlc case, 
heating and cooling by the Compton energy exchange can be calculated <ith the 
coherent (magnetic) scattering cross section (see below). The shape of the 
spectrum should, however, be affected to some degree by incoherent scattering, 
especially at the resonance. The magnetized incoherent scattering is 
numerically very complicated, depending on angle of propagation rather 
sensitively (Nagel 1981b). Under those circumstances where the Compton effect 
represents a small fraction of the total cooling, we expect the spectral 
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distortions to be relatively sma 11 over most of the spectrum. In this paper 
we proceed under this assumption, and compute the continuum spectrum using the 
coherent scattering cross section, although in the energy balance we do 
Include the Compton energy losses and gains. 

For finding the structure of the atmosphere we shall be using the 
expressions for the continuum opacities without the complication of the 
resonance, which is incorporated In the cyclotron line emlsslvlty. Since we 
have a plane parallel atmosphere. It Is reasonable for the purposes of finding 
the atmospheric structure and the spectrum to use the two-stream aproximation 
of the transfer equations, as in Meszaros et al (1980). We use the following 
approximations to the angle averaged, non-resonant opacities in a two 
polarization magnetised plasma: 


S 1 “ "e° 


T + MM + V 2 n e a T AU) 


s 2 = n e a T Kt“ W-) + ( — [l~A(a)>] +V 2 n e a T AU) 
s 12 = S 21 = n e q T ¥ " + ~ "^ +1 /4n e o T A(w) 


a = -_L- [n 2 4tt 2 a 3 
i n e u T L e v th 


(hc/rn ) 2 . -hw/kT 

e — )] 


(7) 


U>‘ 


for w < , 


w-(wm-Au). ) ? 

where u = (u> M /c*> ) 2 , and a(w) * exp [- [ — — =— ) ] 


while for u> > u> H , we take 


; i ■ l h Vt 


0) 
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* V4 Vt 


, _ S 1 r„ 2 4 2 3 (hc/ln e )2 ,1 - <d? /kT 
1 ~ n e a T [ e * » th ** 


•)] • 


( 8 ) 


The Compton scattering opacity, s^, and the bremsstrahlung opacity, a^, are 
all measured in cm -1 , o> is the cyclotron circular frequency 
and v th 5 (irkT/2m) 1//2 . 

We have taken into account the finite width of the resonance with the 
function a(w), which causes the cross sections to approach Oj at w^-Au> l , 
the frequency at which the cyclotron photons are emitted in the lower line 
wing. Cyclotron photons are also emitted in the adjacent frequency bin 
W|_l + Aw^. These opacities are used in performing the two stream radiative 
transfer calculation, layer by layer. 

• • • • 

The bremsstrahlung emissivity is just j’ = a^, where is the Planck 

function for polarization i (ty^of the usual) and a* is the bremsstrahlung 
absorption coefficient when the gaunt factors are set to unity. This 
expression plus the cyclotron line emissivity above is used as the emission 
coefficient in the radiative transfer equations. For the energy balance, we 
use the frequency integrated bremsstrahlung emissivity. 


A b = 2 Q 1 p 2 T ^ 2 (1 - <e^/e^ B >) erg cm " 3 s " 1 (9) 

where we have accounted for inverse effects by including the term <e 1 /e 1 o>, 

V Vd 

the ratio of the actual calculated radiation energy density per unit frequency 
in polarization i to the blackbody energy density e^ B , averaged over all 
frequencies. This approximate form is used in order to have the simpler 
dependence p 2 ^ 2 involving the two atmospheric parameters, this being the 
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frequency Integrated expression of the emissivity. The reason for using this 
simpler expression Is that In the relaxation scheme one has to take numerical 
derivatives of (9) and repeat, to find the equilibrium configuration. The 
values used for are 


Q 1 = 0.6 Q 2 x h -2 [1 - e'*H (1 + x H )] (10) 

Q 2 = 0.4 Q, 


where Q Is the nonmagnetic value 5 x 10 20 (Zeldovich and Shakura 1969) 
and x^ « hv H /kT. These Q - * are the frequency integrated form of the approximate 
(magnetized) bremsstrahlung emlsslvitles = j y NM 0.6 (v/v^) 2 exp ( -hv/kT) 
and j 2 = j 0.4 exp (-hv/kT), j v NM being the nonmagnetic emisslvlty per 
unit frequency. 

The Inverse Compton cooling, in the nonmagnetic case. Involves in the 
nonrelativistic' approximation an expression of the type A - c ay n e U ^ 

(4 kT/m e c 2 ) erg cm" 3 , for electrons at temperature T, U p ^ being a frequency 
integrated photon energy density and cry the coherent Thomson cross section. 

The inverse effect can be taken into account by Including the recoil (eg. 
Kompaneets 1957), so that the average fractional photon energy change per 
scattering in the nonmagnetic case is Av/v - (4kT - hv)/m c 2 . In the strong 
magnetic field case, for hv^ << m e c 2 one can expect that a similar expression 
involving the coherent magnetized scattering cross section will give the 
correct first order Compton heating and cooling. We therefore write the net 
Compton cooling rate in erg cm“ 3 s _1 as 



8<jy^ p Jdv g y (T - hv/4k) 


( 11 ) 
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1 1 

where cr = s /n e from Eqns (7) arid (8). In the approximation where e v Is 
isotropic, one can show from equilibrium considerations that the relative 
value of the terms accounting for direct and Inverse effects should also 
be (T - hv/4k), as In the nonmagnetic case discussed by Ryblckl and Llghtman 
(1979) . 

III. METHOD OF CALCULATION 

We calculate the temperature and density structure of the atmosphere by 
setting up the equations of momentum and energy conservation In the presence 
of the ram pressure and heating by the infalling beam of protons, 

together with the two-stream magnetized 
radiative transfer equations. An initial guess for the atmospheric run of 
temperature and density is made, and the atmosphere is divided into 
homogeneous slabs (20-30) that conform stepwise to these initial guesses. The 
radiative transfer equation is solved layer by layer following the method 
described in the Appendix, and the radiation energy densities for each slab 
are found. We then test each slab for energy balance by comparing the heating 
rate from the accreting protons to the total cooling rate given in Section II 
b: 


* (p ,T,E) = A (p.T.e 1 ), (12) 

where <J> 155 M/A is the proton flux 

and dE/dz = - p dE/dy as a function of (p,T,E) comes from our deceleration 
calculation (stored in tables) . E is an average over the proton energy 
distribution, and a(p, T, e 1 ) contains as a parameter the radiative energy 
densities found by the transfer calculation. We also test for momentum 


balance, 


ar + * *§■ (p,T,E) + (r"+- z -)z” 0 


(13) 


where P = 2n e kT is the thermal gas pressure^ 6 is the gravitational constant, 
and the proton momentum deposition rate dn/dz - - p dn/dy also comes from the 
tables. We have assumed radiation pressure to be unimportant, (this can be 
verified after the solution is found). The proton energy E at the top of the 
atmosphere has the non-relativistic value 


E Q = 134 (-M/M § ) ( 1 0 6 cm/R | j ) MeV, 


(14) 


and the run of £(z) is found by numerically integrating dE(p,T,E)/dz sTab by 
slab until E drops to zero. If E drops to zero inside a particular slab (i.e. 
not at a boundary), then dE/dz, dn/dz are redefined so that the energy and 
momentum with which the proton entered that slab is uniformly distributed over 
the width of that slab, and the stopping length is defined to lie at the 
midpoint of the slab. The last proton leated slab is considered to be the 
bottom of the atmosphere. Each layer, or slab, including the bottom one, is 
assumed to radiate both upwards and downwards. The bottom layer has as 
boundary condition that of total reflection, while the top one has that of no 
incoming radiation from above. These boundary conditions allow us to solve 
the transfer equations layer by layer (see Appendix). They imply also that no 
radiant energy is transported inward into the star (radiated downward from the 
bottom slab), though in reality there may be some loss of energy into the 
stellar surface. 

The equilibrium solution is found iteratively. If momentum and energy 


are not balanced for the density and temperature run of Iteration j, we choose 
a new (j + l)-th gueis of the run of variables by using a relaxation 
technique, a simple description of which is found, for instance, in Clayton 
(1968). One differentiates the. force and energy equations in the variables (p 
and T) to be varied, inserting the energy and momentum deficit found from 
equations (12) and (13) in iteration j, to find Apj +1 , aTj + 1 , which are then 
added to pj, Tj. The proton stopping calculation is repeated to find the new 
heating rate and ram force in each slab. The radiation transfer is also 
repeated each iteration in order to obtain the new for the atmosphere. 
Although the numerically obtained photon energy densities e* also depend on 
T and p, derivatives of cannot be calculated analytically. Differences 
are used to approximate the derivatives of but there is still some degree 
of uncertainty which slows down the convergance. We were, however, able to 
achieve balance in equations (12) and (13), in all cases, to an accuracy of at 
least 10% in each slab and in most cases to much better than 10%. The slabs 
are fixed in space during the iteration so that the az. are constant and the 

J 

yj are recalculated with each change in pj* The top slab is assumed to have a 
density given by the free fall value, 


Pi 


n = 2 x 10 19 (• 

0 10 16 gs“ 1 


) ( IQ 10 cm 2 ^ 


cm -3 , 


(15) 


and this value acts as a boundary condition. The free-fall density is chosen 
as a reasonable upper boundary of the atmosphere since above this point the 
density of the static atmosphere would fall below that of the infalling proton 
beam. This choice does not assign a density to a fixed distance above the 
stellar surface, because the spatial extent of the atmosphere, z Q , is 
determined by the point at which the proton beam stops. Since the protons do 


not necessarily stop in the last slab of our numerical grid, the bottom of the 
atmosphere varies during the iteration. The variations in slab (j + 1) are 
found as a function of those of slab j and the energy and momentum deficits, 
so that the atmosphere balances from top to bottom. The final spectrum Is 
then calculated for the equilibrium atmosphere, using the exact cross 
sections. 

The calculation of the beaming of the radiation requires an angle 
dependent transfer scheme. We use the Integral equation method for a 
magnetized medium (Meszaros and Bonazzola 1981), Knowing the temperature and 
density of each slab, we used the upward directed, angle-dependent intensity 
from each slab as the input for the slab above which then rescatters this beam 
of radiation, as well as itself contributing some emission and absorption. 

The scattering is assumed to be coherent and the bottom slab was assumed to be 
only self-radiating. This method is expedient, in view of the numerical 
magnitude of the task, but it involves the approximation of neglecting the 
backward-scattered radiation. For most of the frequency and magnetic field 
values involved in these calculations, the dominant polarization is 
extraordinary, for which t _< 1. In this case, the approximation of neglecting 
backscattered radiation does not introduce a large error. In addition, a 
property of the magnetized scattering cross sections is that they approximate 
fairly well the case of complete angular redistribution (Nagel 1981a). For 
those frequencies where t > 1 , complete redistribution ensures that scattered 
photons lose Information about their previous directionality, so that the 
inaccuracy incurred is mostly in the total Intensity calculated at the top, 
but not in the calculated angular structure. For this reason, the beam 
structure and pulse profiles are normalized to unit intensity, while the total 
intensity and spectra are calculated by the two-stream approximation, which 


has a full treatment of the boundary conditions. 


IV. RESULTS 

The free parameters determining the structure of the atmosphere are the 

accretion rate ft, the neutron star mass M and radius R^, the surface magnetic 

field strength B and the escape probability of line photons P esc . The 

procedure described in Section IV has been carried out for four accretion 

rates, ft = 10 14 , 1 0 1 5 , 1 0 1 6 , 10 17 g s~ 1 and four values of the magnetic field 

strength B = 2.5xlQi 2 G, 5 x 10 12 G, 10 13 G, 2xl0 13 G,with the other parameters for 

the most part held constant at the values M s 1 M ft , R,,= 1.2 x TO 6 cm., P = 1. 

w n esc 

Figure 2 shows the equilibrium temperature and density profiles for an 

accretion rate of ft = 10 15 g s’* 1 and four different values of magnetic field 

strength. The cyclotron resonance frequency, depending on the magnetic field 

as = 12 B 12 keV, acts as a thermostat which determines the equilibrium 

value of the electron temperature. As a result of the exponential dependence 

of the cyclotron photon production rate on T [cf. eqn. (5)3, cyclotron cooling 

becomes very effective when temperatures in the atmosphere approach 

kT = The equilibrium temperatures are therefore extremely sensitive to 

B and are always less than iiw^. At the bottom of the atmosphere {large y), 

bremsstrahlung and cyclotron line emission are the dominant cooling 

mechanisms, with the relative contributions a function of B. For 
12 

B = 2.5x10 G, cyclotron emission makes up less than 10% of the total 
cooling, whereas for B = 2xl0 13 G v it contributes nearly 90%, with a sharp 
increase in the cyclotron cooling contribution between 5xl0 12 G and 10 13 G . 
Although this result would seem at first to contradict the behavior of the 
cyclotron line cooling expression [cf. eqn. (5)] as a function of B, the 
relative contribution of line cooling actually depends on the ratio of the 
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equilibrium kT to-fta^. As It turns out, equilibrium atmospheric temperatures 
rise with Increasing B faster than the corresponding rise in-hu^, The 

kT/tiw^, therefore Increases with Increasing B, giving a larger line 
cooling contribution. At the top of the atmosphere (small y), Compton cooling 
dominates, due to the dependence on n e rather then n|, as is the case for 
bremsstrahlung and cyclotron cooling. Compton heating makes a significant 
contribution to the energy balance for the higher field strengths, where a 
large number of cyclotron line photons are produced with energies higher than 
that of the atmospheric electrons. When these photons scatter on their way 
out of the atmosphere, they give a portion of their energy back to the 
electrons. 

The density profiles have a roughly linear dependence on y (exponential 
in z) which is characteristic of an isothermal atmosphere in hydrostatic 
equilibrium, with departures duo to the temperature gradients in the 
atmosphere. Scale heights vary from 50 cm. for B = 2.5xl0 12 G to 200 cm. 
for B = 2xl0 13 G. The depth of the atmosphere is the stopping distance of 
the protons, y 0 . The stopping distances decrease with increasing magnetic 
field mainly because the proton energy deposition rates increase with T (cf. 
Fig. 1). Since the optical depth of the atmospheres are directly proportional 
to the stopping distance, 

t - ay o /m p - 6(0/O.5ay) (y o /30 g cm ^), (16) 

atmospheres with higher magnetic fields should be more optically thin. 

This is indeed the case, as is shown in Figure 3, where we have plotted 
the spectrum of radiation emerging from the top of the atmospheres of Figure 
2. The lower field strength photon spectra have a Rayleigh-Jeans behavior up 
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to around 3 keV, which Is the equivalent blackbody temperature, and fall off 
steeply at higher energies, indicating higher optical depths. The atmospheres 
with B > 5xl0 12 G have spectra which begin to resemble optically thin thermal 
bremsstrahlung: power laws of Index ~1 with high energy cutoffs. Figure 3 
shows only the continuum part of the emergent spectrum. Due to our 
approximate treatment of the cyclotron line radiation transfer, we have no 
accurate Information on the appearance of an observable emission (or 
absorption) line in the spectrum. Although the calculation gives the number 
of photons appearing in the line wings, it does not follow the energy 
redistribution of these line photons as they scatter into the adjacent 
continuum. A fraction of the line photons may diffuse down to (at most) 
frequencies Aw ~ k<T e >. We therefore show only the spectrum of the continuum 
photons, keeping in mind that above - .5 this gives only a lower limit on 
the photon number. 

The contribution of the two polarization modes to the continuum spectrum 
is shown in Figure 4, which includes the vacuum polarization effect as do all 
other cases. This effect is due to the presence of virtual e + e - pairs induced 
by the strong magnetic field (MeszSros and Ventura 1978, 1979), which cause 
the polarization eigenmodes to change from circular to linear at the energy 
* 3 keV (n Q /10 22 cm- 3 ) ^ (4x1 0 12 G/B ) . Above this energy, the vacuum 
dominates over the plasma in determining how radiation propagates in the 
medium and causes both modes to be resonant at the cyclotron frequency 
if w v < At low frequencies, (w<u> v ) and at w - where the plasma 

dominates, we can label the modes in the conventional manner with mode 1 being 
the extraordinary and mode 2 being the ordinary mode. In the vacuum dominated 
region, w < to < modes 1 and 2 do not correspond to the extraordinary and 
ordinary modes as usually defined. The convention used here is the same as in 
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Meszaros and Bonazzola (1981). The "vacuum feature" at 4 , found by Venture 
et, al. (1979) in spectra computed for homogeneous atmospheres, does not 
appear in these spectra because of the density gradients, which spread ^o) y 
over a range of frequencies. If, however, each polarization were measured 
separately, one would see different modes predominate above' and below w v , as 
shown in Fig. 4. 

The atmospheric structure and spectrum is much less sensitive to the 
value of the accretion rate. Figure 5 shows the temperature and density 
profiles for four values of ft with B held constant at SxlO 1 '-G. In contrast 
to the large variation of temperature with varying B (cf. Fig. 2), the 
temperatures change only slightly with varying ft. The stopping distances are 
similar for all four cases, with a small increase (y = 20-30 g cm- 2 ) with 
increasing ft . The main effect of increasing ft is to increase the density at 
the bottom of the atmosphere. There is also, therefore, a decrease in the 

scale height, s'ince y Q is roughly constant, to a value of only * 10 cm. 

1 

for ft = 10 17 g s . Although we have neglected the effects of radiation 
pressure in calculating the atmospheric structure, we are able to compute the 
value of the radiation pressure in the equilibrium atmospheres from the 
radiative transfer. We find that it is small relative to the gas pressure 
for ft - 10 16 g s’ 1 , but estimate that it becomes comparable to the gas 
pressure at the top of the atmosphere for ft = 10 g s . Therefore, we may 
not be justified in neglecting radiation pressure for ft = 10 17 g s’ 1 , although 
for the lower accretion rates, the assumption seems to be valid. An accurate 
calculation of the radiation pressure at high ft would require a more detailed 
treatment of the line photon frequency redistribution. 

Figure 6 shows that there is very little difference in the continuum 
spectra emerging from the atmospheres with varying ft, at least in the 
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approximation of coherent scattering. This is not surprising, since It is the 
optical ctepth (proportional to the equilibrium value of y 0 ) which determines 
the shape of the spectrum and y 0 does not vary much with accretion rate. The 
relative contribution of cyclotron line emission to the cooling does vary 
significantly with ft. For ft * 10 14 g s , line cooling makes up 11% of the 
total cooling, whereas for ft « 10* 6 g s" 1 , it represents almost 95% of the 
total. Compton heating makes a significant negative contribution to the total 
for ft - 10 16 g s'* 1 , so that the spectra for the highest accretion rates may 
look different with the Inclusion of incoherent scattering. Since the photons 
which are responsible for most of the Compton heating are line photons, these 
photons would lose energy in the scattering process and enhance the continuum 
below the line. This effect could produce significant differences in the 
spectrum as a function of accretion rate, at the higher ft values. 

We have so far presented results for cases where P * 1, meaning that 
line photons have a 100% probability of escaping without aosorptlon. Cases 
where P s 0.1 were also Investigated and were not found to produce 

G 

significantly different results from the P aen - 1 cases for the same values 
of B and ft. The reason this is so lies in the approximate expf-tfu^/kT) 
dependence of the line photon production rate. It is possible for a small 
Increase in T to counteract an order of magnitude difference in P esc , such 
that the line cooling remains the same (ie. ten times more line photons are 
produced but only one tenth of them escape). The increase in T necessary to 
do this is small enough to produce little observable difference in the T 
profile or the spectrum. 

We have plotted in Figure 7 the pulse shape, normalized in each case to 
unit intensity, for different values of B, frequency and aspect angles 

These latter are the angle between the magnetic and rotation axis, 
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and that between the line of sight and the rotation axis. The pulse shapes 
are symmetric to an Interchange of <j>-j and <j> 2 . The pulse shapes, produced by 
the radiation passing through atmospheres of decreasing density and slowly 
Increasing temperature, can be compared with those coming from homogeneous 
atmospheres. The present pulses resemble qualitatively those coming from a 
uniform, self-radiating slab, with no background Illumination, discussed by 
Meszaros and Bonazzola (1981). Of course in those previous calculations, the 
choice of temperature and density was arbitrary. One can see, however, that 
different temperature gradients produce different pulse shapes, so that a 
uniform atmosphere at a density and temperature of the bottom slab would still 
be only a rough approximation. 

The pencil beam pulses from these Coulomb heated atmospheres are 
characterized over a significant range of aspect angles by single pulse shapes 
at high frequencies and notched or multiple pulse shapes at low frequencies. 
The energy at which the pulses undergo this metamorphosis Is -tlw •'"tfoi^/4. The 
pulses therefore become multiple at higher energies for higher field 
strengths, where Is larger (cf. Fig. 7). The notches occur at phase zero, 
when the observer Is looking down the magnetic pole. At this angle, the 
magnetized cross sections undergo a minimum at frequencies below If the 

optical depth of the atmosphere Is not too large, the outgolnq intensity will 
be directly proportional to the emlsslvlty which has the same angle and 
frequency dependence as the cross sections. The Intensity will thus have a 
minimum at phase 0. If the atmosphere Is optically thick, the Intensity is 
Inversely proportional to the cross section and the notches tend to disappear 
(Nagel 1981a), The notches In the higher magnetic field cases are therefore 
more pronounced as a result of the lower optical depth. 

The spectral Index is a function of phase, since the Intrinsic beam shape 


is both angle and frequency dependent. As previously found in homogeneous 
slab models (M6szaros and Bonazzola 1981, Nagel 1981a), the spectral Indices 
for most aspect angles are hardest at midpulse. In Figure 8, we have plotted 
the spectral index as a function of pulse phase for the pulses In Figure 7. 
The spectra for B * 5xl0 12 G show a hardening toward phase 0 (midpulse) for 
the energy range 10-18 keV, where the pulses are developing notches, but none 
in the energy range 3-10 keV, where the pulse shapes are stable. Similarly, 
the specta for B * 10 13 G show a hardening toward midpulse only in the lower 
energy range where the pulse notches are deepening. 


V, DISCUSSION 


We have investigated the structure of accreting atmospheres where the 
main deceleration mechanism of the infalling protons is distant Coulomb 
encounters with atmospheric electrons. We assume that a colli sionl ess shock 
does not arise, although a priori one can only make plausibility arguments 
(e.g. Kirk and Trumper 1982). The alternative approach of calculating an 
accreting magnetized atmosphere when a collisionvess shock does occur has been 
explored by Langer and Rappaport (1982), while the effects of radiation 
pressure Including approximate magnetic cross sections was explored by Wang 
and Frank (1981). Earlier work by Basko and Sunyaev (1976) assumed a magnetic 
field to determi/!'- the geometry of the accretion column, but used nonmagnetic 
cross sections. These atmospheres with a collisionless shock or with strong 
radiation pressure (which produces a quasi-shock) differ radically from what 
we have calculated here, since the shock stands off from the stellar surface 
at a distance az which is comparable to or larger than the transverse 
dimension of the accretion column, The optical depth across B can be less 
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than that along B, so that radiation escapes from the sides of the column, 
forming fan beam pulses. The Coulomb heated atmospheres, on the other hand, 
have Az much less than the polar cap radius, which allows us to use the 
convenient plane parallel approximation. 

We find that Coulomb heated atmospheres are indeed thin slabs with scale 
heights of only - 10-200 cm., much smaller than the polar cap 
radii R p - 1 0 4 -1 0 s cm. determined from the accretion flow. Very little 
radiation would be expected to escape from the sides of the slab, since the 
optical depths are ~ 1000 times greater in this direction than along the 
magnetic pole. Densities are quite high (10 2z - 10 24 cm" ) at the base of 
the atmosphere, with the highest densities occurring for low field strengths 
and high accretion rates. 

-2 

The equilibrium stopping distances are In the range y Q « 20 - 50 g cm 
with the longer stopping distances occurring for low field strengths and low 

o 

accretion rates. In the cases where y Q - 50 g cm , nuclear collisions, which 

were included in the energy deposition rate, play a significant role in 

slowing down the protons in the accretion flow. If the Coulomb energy 

deposition rates used in the present calculation are overestimates, then 

nuclear collisions could be expected to be important for the higher field 

o 

strength cases also, the values of y 0 would be higher, - 50 g cm , and the 
continuum spectra would be more optically thick. The spectra for the low 
field strengths would not change much since these cases have y Q close 
to 50 g cm- 2 even with the present Coulomb rates. 

Cyclotron line cooling is very important in determining the equilibrium 
temperature of the atmosphere. Due to the - exp {-liu^/kT) dependence in the 
production rate of line photons by colli si onal excitations, the equilibrium 
temperatures are extremely sensitive to the magnetic field strength. As the 
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field increases, line cooling contributes more and more to the total cooling, 
but becomes less effective in keeping the temperature low. As moves up 
into the > 200 keV range, as is the case for B - 2xl0 13 G, equilibrium 
temperatures approach 10 9 °K, where pair production (either photon-photon or 
photon-magnetic field) could become an important cooling mechanism. It is 
therefore unlikely that in these equilibrium atmospheres, temperatures could 
greatly exceed 10 9 °K. 

We have neglected the excitation of electrons to Landau levels higher 
than the first excited state. This could become important at low field 
strengths, if the energy and spacing of the electron states are small compared 
to the equilibrium temperature in the atmosphere. It is even possible in very 
low fields for the first few cyclotron harmonics to become optically thick, so that 
the cooling contribution would come from a higher transition (the lowest 
optically thin transition). However, for the range of field strengths we have 
considered, the temperatures are all well below the energy of the first 
excited state. Furthermore, radiative ^excitation is so rapid compared to 
colli si onal processes that very few electrons would be expected to be excited 
from the first to the second excited state. In addition, bremsstrahlung is so 
effective for low field strengths that even for fields below 2.5xl0 12 G, there 
may be no need to consider cylcotron cooling from any but the lowest excited 
state. 

We have also neglected the changes in energy of the photons as they 
Compton scatter, though we do take into account the Comptonization of the 
electrons in the cooling part of the energy balance. Although the change in 
energy of a photon is small in each scattering event, each photon may undergo 
many scatterings in the atmosphere before it escapes. One indicator of the 
importance of these incoherence effects is the contribution of Compton 
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scattering to the total cooling, since the energy gained or lost by the 
electrons must equal that of the photons. For the equilibrium atmospheres we 
have presented, Compton scattering makes up about 30% of the cooling or 
heating, except for ft - 10 16 g s" 1 where the percentage Is somewhat higher. 
However, most of this contribution is heating by line photons and not 
scattering of continuum photons. The cooling by lower energy continuum 
photons at the top of the atmosphere represents only a few percent of the 
total cooling so the continuum spectrum below - .5 is unlikely to be 
seriously affected. The main consequences of neglecting incoherence effects, 
therefore, are poor information on the shape of the line, and an underestimate 
of the continuum level Immediately below the line. We may also have somewhat 
overestimated Compton heating by line photons in preventing these photons from 
losing energy as they scatter. As they blend into the continuum they should 
become less effective in heating electrons. 

A number of observed X-ray pulsars have hard power law spectra with 
turnovers at around 10-20 keV (Rose et al. 1979, Pravdo et al. 1979, Becker et 
al. 1977, Pravdo et al. 1978). In several sources, such as 0A01653-40 and 
4U1 145-61, the power law spectra show evidence only for a slight steepening 
above 20 keV, indicating the presence of higher temperatures. Many of these 
pulsars show* broad line features in their spectra at energies between 6.6 and 
7 keV which have been identified as Fe emission. Only two pulsars, however, 
show clear evidence for cyclotron line features in their spectrum: Her X-l 

and 4U0115+63 (Voges et al. 1982, White, Swank and Holt 1983). Our 
calculations suggest that some of the pulsars having hard power law spectra 
without cyclotron features may have higher magnetic fields. A spectrum with 
index -1 extending to 30-50 keV before steepening requires B > 5xl0 12 G; the 
cyclotron frequency would be at > 60 keV, above the sensitivity of most 
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detectors. The spectra for B - 5xl0 12 G fall off steeply in the 10-30 keV 
range, although the downscattering of resonance photons could change the shape 
of the spectrum In this region. 

The best observational evidence fov 1 slab atmospheres which produce pencil 
beam pulses comes from the observed frequency dependence of pulse shapes. 

Many X-ray pulsors show transitions from single pulses at high frequency to 
multiple pulses at low frequencies, without phase shifts (Pravdo et. al. 1977, 
White et. al. 1983), as well as a spectral hardening towards midpulse. Since 
our calculations predict this transition to occur at - 1/4 - 1/3^, an 
independent observational estimate of B is possible. This is a higher 
transition frequency than predicted by Meszaros and Bonazzola (1981) for semi- 
infinite homogeneous atmospheres. The exact fraction of the cyclotron 
frequency below which the pulse splits depends on the viewing angle. 

Since we now have spectral, pulse shape and phase-spectroscopic 
predictions from one of the two main X-ray pulsar models (for - IQ 37 erg. s— 1 ) , 
the time seems ripe for a more detailed comparison with the data. From the 
preliminary comparison we have made, there seem to be several overall 
similarities between some X-ray pulsars and our models, which is 
encouraging. One should bear in mind, however, that there are departures from 
the "average" observed behavior, in addition to the fact that the range of 
applicability of our model, L < 10 37 erg s— 1 , reduces the total available 

A 

sample. 

For the alternative shock models, a spectrum for a self-consistent column 
atmosphere has been given by Langer and Rappaport (1982), which may also fit 
some observed X-ray pulsar spectra. While pulse shape and spectral index vs. 
phase predictions have not yet been made for the shock models, previous work 
on non self-consistent column atmc .spheres suggests that the resulting fan-beam 
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pulses may be rather broad (Yahel 1980, Pravdo and Bussard 1981, Nagel 
1981a). This may be Inconsistent with many pulsars, but determinations of the 
pulse shapes for self-consistent shock atmospheres are needed before one can 
decide on this Issue. 

As the volume of theoretic*. 1 predictions from the shock and Coulomb 
models grows along with observational data on X-ray pulsars, statistical tests 
of consistency may be very fruitful. In this respect, the organization of 
available observational material on spectra, pulse shape and pulse phase 
spectroscopy on as many X-ray pulsars as possible, as provided by White et al. 
1983, is of great value. We may then be able to decide, from consistency 
tests, in favor of one model or the other. Detailed comparison with model 
predictions could furthermore yield extremely useful information on physical 
parameters in the accreting atmospheres, such as magnetic field strength, 
temperature and density. 
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APPENDIX 

RADIATIVE TRANSFER IN AN INHOMOGENEOUS ATMOSPHERE 

Meszaros, Nagel and Ventura (1980) have given solutions of the transfer 
equations for polarized radiation in a strongly magnetized, self radiating 
homogeneous slab. We use these results to numerically solve the one- 
dimensional transfer problem in the two-stream approximation for an 
inhomogeneous plane parallel atmosphere. As discussed In Section V, the 
Coulomb heated atmospheres have scale heights which are much smaller than 
their cross-sectional diameters and so both the plane parallel and the two 
stream approximations should be quite good. 

We denote the upgoing and downgoing intensities in the orthogonal 
polarization modes 1 and 2 by J 2 , , K 2 respectively. The slab has a 

width z 0 , with z = 0 at the stellar surface, and is divided into N layers, 
each having a constant temperature, density, and fixed width Az n< The 
solutions to the transfer equations which are given by Meszaros et al (1980) 
give the intensities in each layer. The solutions in layer n at each 
frequency can be written in the matrix form: 

t n (z) = t n + A n (z) t n (Al) 


where 
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The £ are undetermined coefficients, the 4 x 4 matrices A p (z) involve 
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combinations of scattering and absorption coefficients multiplied by 
exponentials In z, and | is the free-free and the cyclotron line emission in 
the layer. The solutions In each layer Involve 4 unknown coefficients; and N 
layers minus 4 boundary conditions gives a total of 4N-4 unknowns. We solve 
for these 4{N-1) coefficients simultaneously by matching the Intensities at 
each of the (N-l) boundaries. The boundary conditions on the intensities at 
z=z n can then be written, 

£ n -l + A n-1 <z n* ^n-1 = ^n + A n (z n } ^n (A2) 


As discussed In Harding and Tademaru (1981), simultaneous equations of this 
form can be solved without inverting N X N matrices, by writing the 
coefficients In layer N in terms of those in layer 1 as 


Afj ( Z N ) C N = % Aj (z 2 ) + t H 


(A3) 


where 


+ N N-l 

X N = ii=2 nfln A m (z m+l 



^n-1"V 


and 


M N ■ & A n <Vl> A n’ < z n>- 


We choose the boundary conditions to be 


1 _ 

1,2 " l ,2* 


(A4) 


That Is, no radiation Is Incident on the slab from above, and the downgoing 
radiation In the bottom slab Is totally reflected. With these boundary 
conditions, the problem of finding the Intensities emerging from the top and 
bottom of the slab has been reduced to solving the set of four simultaneous 
equations (A3). The Intensities In all other layers can then easily be 
determined from the bottom layer up using equation (A2). 
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FIGURE CAPTIONS 


Figure 1: 


Figure 2: 


Figure 3: 


Figure 4: 


Figure 5: 


Figure 6: 


Proton energy loss rate as a function of atmospheric electron 
temperature for different proton energies In MeV. 

Atmospheric temperature and density profiles for different 

magnetic field strengths In units of 10 12 G, plotted against y * 
z 

- / njz)m dz. Infall Is from the left, 
oo e P 

Photon number spectra for the atmospheres of Fig. 2, labeled with 
different values of the magnetic field strength. Only the 
continuum spectra are shown, although additional photons are also 
produced at the cyclotron frequency, which has the values 30, 60, 
120 and 240 keV for these cases. 

Photon number spectrum for the case B = 5xl0 12 G, showing the 
spectra of individual polarization modes (1 = extraordinary, 2 = 
ordinary). 

Atmospheric temperature and density profiles for different 
accretion rates and constant field strength B = 5x10 12 G. 

Photon number spectra for the atmospheres of Fig. 5, labeled with 
different values of the accretion rate in g s“*. The cyclotron 
frequency for all cases is 60 keV. 


FI gun 7: 


Normalized pulse shapes for aspect angles (60°, 45°) 
ft * 10 15 g s” 1 , P esc * 1, and four different values of magnetic 
field strength, each at three photon energies: 18, 10, and 3 keV 
from top to bottom, 

Figure 8: Energy spectral Index a, defined by dL s KE a dE, as a function of 

pulse phase, in two different energy ranges: 3-10 keV (solid 
lines) and 10-18 keV ( dashed lines), for aspect angles: 
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